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ABSTRACT 

In the braneworld scenario the four dimensional effective Einstein equation has ex- 
tra source terms, which arise from the embedding of the 3-brane in the bulk. These 
non-local effects, generated by the free gravitational field of the bulk, may provide an 
explanation for the dynamics of the neutral hydrogen clouds at large distances from 
the galactic center, which is usually explained by postulating the existence of the dark 
matter. In the present paper we consider the asymptotic behavior of the galactic ro- 
tation curves in the brane world models, and we compare the theoretical results with 
observations of both High Surface Brightness and Low Surface Brightness galaxies. 
For the chosen sample of galaxies we determine first the baryonic parameters by fit- 
ting the photometric data to the adopted galaxy model; then we test the hypothesis 
of the Weyl fluid acting as dark matter on the chosen sample of spiral galaxies by 
fitting the tangential velocity equation of the combined baryonic- Weyl model to the 
rotation curves. We give an analytical expression for the rotational velocity of a test 
particle on a stable circular orbit in the exterior region to a galaxy, with Weyl fluid 
contributions included. The model parameter ranges for which the test provides 
agreement (within la confidence level) with observations on the velocity fields of the 
chosen galaxy sample are then determined. There is a good agreement between the 
theoretical predictions and observations, showing that extra-dimensional models can 
be effectively used as a viable alternative to the standard dark matter paradigm. 

Key words: cosmology: dark matter - galaxies: haloes - gravitation: relativistic 
processes. 



1 INTRODUCTION 

Gravitational effects which require more matter than what 
is visible can be explained in terms of a mysterious dark 
matter, the nature of which remains a long-standing problem 
in modern astrophysics. Two important observational issues, 
the behavior of the galactic rotation curves and the mass 
discrepancy in clusters of galaxies, led to the necessity of 
considering the existence of dark matter both at galactic 
and extra-galactic scales. 

T he rotation c ur ves o f spiral ealax- 
ies (iBinnev fc Tremaiii^ Il987l : iPersic et al.1 Il996l : 
iBoriello fc Saluccil bOOlf T are among the best evidences 
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showing the problems Newtonian gravity and/or standard 
general relativity have to face on the galactic/intergalactic 
scale. In these galaxies neutral hydrogen clouds are observed 
at large distances from the center, much beyond the extent 
of the luminous matter. Since these clouds move in circular 
orbits with velocity vtg{r), the orbits are maintained by the 
balance between the centrifugal acceleration Vtg/r and the 
gravitational attraction GM(r)/r'^ of the total mass M{r) 
contained within the orbit. This allows the expression of 
the mass profile of the galaxy in the form M{r) = rvlgjG. 

Observations show that the rotational velocities in- 
crease near the center of the galaxy, in agreement with 
the theory, but then rem ain nearly co nstant at a value of 
vtgoo ~ 200 - 300 km/s (|Binnev fc Trem aine 1987), which 
leads to a mass profile M(r) = rvtgoo/G. Consequently, the 
mass within a distance r from the center of the galaxy in- 
creases linearly with r, even at large distances where very 
little luminous matter has been detected. 
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The second astrophysical evidence for dark matter 
comes from the study of the clusters of galaxies. The to- 
tal mass of a cluster can be estimated in two way s. Know- 
ing th e motions of its member galaxies ,Giovanelli et al.l 
l|l994 ). the virial theorem gives one estimate, Mvt, say. 
The second is obtained by separately estimating the mass 
of the individual members, and summing up these masses, 
to give a total baryonic mass Mb- Almost without excep- 
tion it is found that Mvt is considerably larger than A'Ib, 
Mv t > Mb, typical values of Mvt /Mb being about 20- 
30 l|Binnev fc Tremaine|[l987l ). 

This behavior of the galactic rotation curves and of 
the virial mass of galaxy clusters is usually explained by 
postulating the existence of some dark (invisible) matter, 
distributed in a spherical halo around the galaxies. The 
dark matter is assumed to be a cold, pressure-less medium. 
There are many possible candidates for dark matter, the 
most popular ones being the weakly interacting massive 
particles (WIMP) (for a revie w of the particle physics as - 
pects of the dark matter see lOverduin fc WessonI l|2004h . 
Their interaction cross sections with normal baryonic mat- 
ter, while extremely small, are expected to be non-zero, 
and we may expect to detect them directly. Models based 
on right-handed (sterile) neutrinos (|Biermann fc Kusenkol 
l2006l : fMunvaneza fc Biermannll2006h as a form of warm dark 
matter were also advanced. 

It has also been suggested that the dark matter in the 
Universe might be composed of super-heavy particles, with 
mass ^ 10^" GeV. But observational results show that the 
dark matter can be composed of super-heavy particles only 
if these interact weak ly with normal matter, or if th eir mass 
is above 10^^ GeV l| Albuquerque fc Baudid l2003h . Scalar 
fields or other long range coherent fields coupled to gravity 
have also intensively been used to model galactic dark mat- 
ter ( Fuchs fc Miclkc 2004; Hernandez ct al. 2004; Gia nni^ 
20051 : [Guzman fc Urena-Loped l2006l : iBernal fc GuzmanI 



20061 ; iBriscese .20111 : iHarkdbOllI ) 



However, up to now no non-gravitational evidence for 
dark matter has been found, and no direct evidence or an- 
nihilation radiation from it has been observed yet. 

Therefore, it seems that the possibility that Einstein's 
(and the Newtonian) theory of gravity breaks down at the 
scale of galaxies cannot be excluded a priori. Several theo- 
retical models, based on a modification of Newton's law or of 
general relativity, have been propos ed to expl ain the be hav- 
ior o f the galactic rotation cur ves llMilgrom| [r983: Sander^ 



1984 



2004 



Moffat fc Sokolovl Il996l: |Mannheii ij| [1997I: iRoberti 
Bohmer fc Harkd l2007al: Bohincr fc Harko I l2007b| : 



Bcrtolami ct al. 2007; Boh mer et al.ll2008inBekensteinll2004 ; 
Brownstein fc Moffat .2006 ). 

The idea of embedding our Universe in a higher 
dimensional space has attracted a considerable inter- 
est recent ly, due to the proposal by Randall and 
Sundrum (|Randall fc SundrumI Il999al lbh that our four- 
dimensional (4D) space-time is a three-brane, embedded in a 
5D space-time (the bulk) . According to the brane world sce- 
nario, the physical fields (electromagnetic, Yang-Mills etc.) 
in our 4D Universe are confined to the three brane. Only 
gravity can freely propagate in both the brane and bulk 
space-time, with the gravitational self-couplings not signifi- 
cantly modified. Even if the fifth dimension is uncompacti- 
fied, standard 4D gravity is reproduced on the brane in the 



appropriate limit. Hence this model allows the presence of 
large, or even infinite non-compact extra dimensions. Our 
brane is identified to a domain wall in a 5D anti-de Sitter 
space-time. For a revie w of the dynamics and geom etry of 
brane universes, see e.g. lMaartens fc Kovamal (|201G| ). In the 
brane world scenario, the fundamental scale of gravity is 
not the Planck scale, but another scale which may be at the 
TeV level. The gravitons propagating through the bulk space 
give rise to a Kaluza-Klein tower of massive gravitons on 
the brane. These gravitons couple to the energy-momentum 
term of the standard model fields, and could be produced 
under the appropriate circumstances as real or virtual par- 
ticles. 

Due to the correction terms coming from the extra 
dimensions, significant deviations from the standard Ein- 
stein theo2;_occur in brane w orld mode ls at very hi gh ener- 
gies (ISasaki et al., 2 000: Shi rornizu et al. 2000: Macd a~et al.l 
l2003l : lGergelvll2003l , |2008| . |2009| ) . Gravitv is largely modified 
at the electro-weak scale of about 1 TeV. The cosmological 
and astrophysical implications of the brane world theories 
have been extensively investigated in the physical literature. 

The static vacuum gravitational field equations on the 
brane depend on the generally unknown Weyl stresses, which 
can be expressed in terms of two functions, called the dark 
radiation U and the dark pressure P terms (the projections 
of the We yl curvature of the bu l k, generating non-local brane 



or tne we yl curvature oi tne du i k, generatmg non-iocai prane 
stresses) IIDadhich et al. ll2000l:lGermani fc Maartendl 200 j; 



is) 

Gergelvll2006l : iBohmer fc Harkdl2007d : iMaartens fc Kovamal 



20101 ). Generally, the vacuum field equations on the brane 



can be reduced to a system of two ordinary differential equa- 
tions, which describe all the geometric properties of the vac- 
uum a s functions of the d ark pressure and dark radiation 
terms (|Harko fc Mak|[2003t ). In order to close the system of 
vacuum field equations on the brane a functional relation 
between these two quantities is necessary. 

The full 5-dimensional Einstein equations have been 
solved numerically for stat ic, spherically sy mmetric matter 
localized on the brane in IWisema n| (12002), yielding regu- 
lar geometries in the bulk with axial symmetry. The same 
data that specifies stars in 4-dimensional gravity, uniquely 
constructs a 5-dimensional solution. An upper mass limit 
is observed for these small stars, and the geometry shows 
no global pathologies. The intrinsic geometry of large stars, 
with radius several times the AdS length, is described by 
four-dimensiona l General Relativity. The results obtained in 
|w isema 3 (|2002D show that Randall- Sundrum gravity, with 
localized brane matter, reproduces relativistic astrophysical 
solutions, such as neutron stars and massive black holes, in 
a way which is consistent with observations. 

Several classes of spherically symmetric so- 
lutions of the static gravitational field equations 
in t he vacuum o n the bran e have been o btaine d 
Harko fc Ma^ (|2003|): iMak fc Harkd h004]: 

iGergelvl 
a possi- 



( 20071 ): iHorvath. Gergely fc HobTIil dlOOsI). As 



ble physical application of these solutions the behavior 
of the angular velocity vtg of the test pa rticles in sta- 
ble cir cular orbits ha s been considered (iMak fc Harkd 



20041: iHarko fc Chend 
Rahaman et al.l bOOSj ) 



20061 : iBohmer fc Harkd 120073 : 
The general form of the solution, 
together with two constants of integration, uniquely deter- 
mines the rotational velocity of the particle. In the limit 
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of large radial distances, and for a particular set of values 
of the integration constants the angular velocity tends to a 
constant value. This behavior is typical for massive particles 
(hyd rogen clouds) outside galaxies (|Binnev fc Tremaind 
1 19871 ). and is usually explained by postulating the existence 
of the dark matter. 

Thus, the rotational galactic curves can be naturally 
explained in brane wo rld models, without introducing any 
addit ional hypothesis jMak fc Harkdl2004l : iHarko fc Chend 
I2OO6I : iBohmer fc Harkdl2007d '). The galaxy is embedded in 
a modified, spherically symmetric geometry, generated by 
the non-zero contribution of the Weyl tensor from the bulk. 
The extra terms, which can be described in terms of the 
dark radiation term U and the dark pressure term P, act as a 
"matter" distribution outside the galaxy. The particles mov- 
ing in this geometry feel the gravitational effects of U, which 
can be expressed in terms of an equivalent mass (the dark 
mass) Mu ■ The dark mass is linearly increasing with the dis- 
tance, and proportion al to the baryonic mass of the galaxy, 
Mu{ r ) K. Mnir/m) jMak fc Harkol |2004| : iHarko fc Chen3 
I2OO6I : iBohmer fc Harkdl2007d '). The exact galactic metric, 



the dark radiation and the dark pressure in the flat rotation 
cu rves region in th e bran e world scenario has been obtained 
in lHarko fc Chend l|2006l ). 

Similar interpretations of the dark matt er as bulk effect s 
have been also considered in iPal et all (|2005l ') and|Pai (|2008D , 
where it was also shown that it is possible to model the X- 
ray profiles of clusters of galaxies, without the need for dark 
matter. 

Hence, a first possible approach to the study of the vac- 
uum brane consists in adopting an explicit equation of state 
for the dark pressure as a function of the dark radiation. 
A second method consists in closing the field equations on 
the brane by imposing the condition of the constancy of the 
rotational velocity curves for particles in stable orbits. 

It is the purpose of the present paper to consider the 
general behavior of the vacuum gravitational field equations 
in the brane world model in the region of constant tangential 
rotational velocity of test particles in stable circular orbits, 
and to compare the predictions of our model with the exist- 
ing observational data. Physically, this situation is charac- 
teristic for part icles gravitating in circular orbits around the 
galactic center (|Binnev fc Tremaindll987h . As a first step in 
our study we derive, under the assumption of spherical sym- 
metry, the basic equations describing the structure of the 
vacuum on the brane, and the equation giving the behavior 
of the tangential velocity of the test particles as functions of 
the dark radiation and of the dark pressure. 

To obtain the tangential velocity of test particles we 
analyze the geodesies for a 4D motion by reducing the prob- 
lem to the motion of a test particle in an effective potential, 
also containing the effects of the bulk. By assuming that the 
brane is a fixed point of the bulk, we derive the equation 
giving the tangential velocity as a function of the gtt com- 
ponent of the metric tensor only. In the region of constant 
tangential velocities, the general solutions of the gravita- 
tional equations can be obtained in an exact analytic form. 

To test the viability of the model, and in order to apply 
it to realistic systems, we compare the theoretical rotation 
curves with a sample of 18 galaxies that includes both low 
and high surface brightness galaxies, with measured rota- 
tion curves extending in the dark matter dominated region. 



From this comparison we obtain the numerical values of the 
parameters characterizing the equation of state of the Weyl 
fiuid, as well as the parameters giving the behavior of the 
rota tion curves i n the b rane world models. 

iPersic et al.1 l|l996l ) found that spiral rotation curves 
show that the rotation velocity at any radius depends only 
on the luminosity. This implies that the galactic rotation 
curves can be modeled by a Universal Rotation Curve. These 
curves show a slightly declining velocity in the galactic halo. 

The present paper is organized as follows. The field 
equations in the brane world models, as well as the tan- 
gential velocity of test particles in stable circular orbits are 
presented in Section 2. The structure of the vacuum in the 
brane world models is analyzed in Section 3. The explicit 
form of the tangential velocity is obtained, for a specific 
form of the equation of state of the Weyl fiuid, in Section 
4. A general theoretical comparison of the tangential veloc- 
ity relation with the observational constraints is performed 
in Section 5. In Section 6, by fitting the theoretical predic- 
tions of the model with the observational data for a sample 
of 9 HSB galaxies, the numerical values of the parameters 
describing the Weyl fiuid are determined. We use the same 
method for the LSB galaxies in Section 7. We discuss and 
conclude our results in Section 8. 



2 FIELD EQUATIONS AND TANGENTIAL 
VELOCITY OF TEST PARTICLES 

In the present Section we present the field equations for 
static, spherically symmetric vacuum branes, and obtain the 
velocity of the test particles in stable circular orbits around 
the galactic center. 



2.1 The field equations in the brane world models 

We start by considering a five dimensional (5D) space-time 
(the bulk), with a single four- dimensional (4d) brane, on 
which ordinary matter is confined, only gravity can probe 
the extra dimensions. The 4D brane world (}^^M,gab) is lo- 
cated at a hypersurface (^B {X^'j = O) in the 5d bulk space- 
time C"^' A'/, (/as), of which coordinates are described by 
X^,A — 0,1,. ..,4. The 4D coordinates on the brane are 
x",a — 0,1,2,3. All tensors Tab and vectors Va satisfy 

T ij = T ahSiSj and Va = VuSa at the hypersurface. 

The action of the system is given by (ISasaki et aLl lioOol) 



S — Sbulk + Sbra 



where 

Sbulk - 

and 



(S) A/ 



(5) J 



2kl 



(1) 



(2) 



-r^K + Lbrane {Qac, 1p) + ^b 



a X, 



(3) 



Sbrane = / V^^^ >< 

where fef = 87rG5 is the 5D gravitational constant, R 
and ^-^^ Lm are the 5D scalar curvature and the matter 
Lagrangian in the bulk, representing non-standard model 
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fields, Ltrane (flat,''/') IS the 4D Lagrangian, which is given 
by a generic fianctional of the brane metric gat and of the 
matter fields ip, is the trace of the extrinsic curvature 
on either side of the brane, and A5 and A;, (the constant 
brane tension) are the negative vacuum energy densities in 
the bulk and on the brane, respectively. 

The Einstein field equations in the bulk can be obtained 
as (|Sasaki et al.ll2000l ) 



(5) 



Gij = kl^'^Ti 



where 



(5) 



Tij = -2- 



in -L (5)„„(5) 



(4) 



(5) 



is the energy- momentum tensor of bulk matter fields. The 
energy-momentum tensor of the matter localized on the 
brane, is defined by 

ane 



T, 



H~ Qab-^brane ■ 



(6) 



(The minus sign in the above definitions appears because 
the variation is done with respect to fl"'' rather then gab, 
and there are two terms because the Lagrangian is varied, 
rather then the Lagrangian density.) 

The delta function 5 (B) denotes the localization of 
brane contribution. In the 5D space-time a brane is a fixed 
point of the Z2 symmetry. The basic equations on the brane 
are obtained by projections onto the brane world. The in- 
duced 4D metric is gu = ^^^gij — ninj, where ni is the 
space-like unit vector field normal to the brane hypersur- 
face ''''M. In the following we assume '^'l/m = 0. 

Assuming a metric of the form ds^ = {ninj + 
gij)dx^ dx'' , with nidx' = dx the unit normal to the 
X = constant hypersurfaces and gij the induced metric 
on X = constant hypersurfaces, the effective 4d gravita- 
tional equation on the brane (the effective Einstein equa- 
tion), takes the form (Sasaki et al..,200Q '): 

Gab = k\Tab + kfSab — E ab, (7) 

where Sab is the local quadratic energy-momentum correc- 
tion 



3T'"'Tcd - T^) 



(8) 



Sab = — riai) — -Ta'^Tijc + 24^"'' 

and E ab is the non-local effect from the free bulk gravita- 
tional field, the transmitted projection of the bulk Weyl ten- 
sor CiAjB, E ij = CiAjBiT-^n^ , with the property E jj 
E abSjSj as X — > 0. We have also denoted fc| = StvG, with 



G the usual 4D gravitational constant. In the limit 







we recover standard general relativity (jSasaki et al.ll200ol l. 

The Einstein equation in the bulk and the Co- 
dazzi equation also imply the conservation of the energy- 
momentum tensor of the matter on the brane, VcTa'^ = 0, 
where Va denotes the brane covariant derivative. Moreover, 
from the contracted Bianchi identities on the brane it fol- 
lows that the projected Weyl tensor obeys the constraint 
VciJa" — kf^cSa"- Thus the 4-di vergence of Eab is con - 
strained by the matter on the brane l|Shiromizu et al.ll2000D . 

The generic traceless Eab can be decomposed ir- 
red ucibly with respect to a chosen 4-velocity field u"' 
as l|Maartens fc KovamalbOloD 



E ab 



u 



(uaUb + ^hab^ 



+ Pab + 2Q 



(a^b) 



(9) 



where the induced metric hab = flat + UaUb projects orthog- 
onal to u° , the "dark radiation" term U = —k^'^ E abu°'u!' 



is a scalar, Qa 
P u — —k'"- 



k. 



^^'^ ha E cdu"^ is a spatial vector, finally 
^ab = — 'C4 [h'^a^b) ~ l^^at'i'^'^] -E cd IS a Spatial, Symmetric 
and trace- free tensor. 

In the following we neglect the effect of the cosmologi- 
cal constant on the geometry and dynamics of the galactic 
particles. In the case of the vacuum state we have p — p — 0, 
Tfj^ = 0, and consequently Sab ~ 0. Therefore the field equa- 
tion describing a static brane takes the form 



Rab — ~ E ab 



(10) 



with the trace R of the Ricci tensor Rab satisfying the con- 
dition R = R'^ = 0. 

In the vacuum case E ab also satisfies the constraint 
VcEa'^ ~ 0. In an inertial frame at any point on the brane 
we have it" — 5q and hab = diag(0, 1, 1, 1) . In a static 
vacuu m Qg = and the const raint for E ab takes the 
form (|Germani fc Maarten^l 20011 ) 



^DaU + ^UAa + D'Pab + A" Pab = 0, 



(11) 



where Aa — u VbUa is the 4-acceleration and Da denotes 
the covariant derivative associated to the metric hab- In 
the static spherically symmetric case we may chose Aa — 
A(r)ra and Pab = -P(r) (r^r;, - ^hab), where A{r) and P{r) 
(the "dark pressure" although the name dark anisotropic 
stress might be more appropriate) are some scalar func- 
tion s of the radial dista nce r, and ra is a unit radial vec- 
tor (|Dadhich et al.ll2000h . Thus the expression (|9]) simplifies 
to 



E, 



-ki 



UUaUb + 



U-P 



hab + PVaTb 



(12) 



2.2 The motion of particles in stable circular 
orbits on the brane 

In brane world models test particles are confined to the 
brane. Mathematically, this means that the equations gov- 
erning the motion are the standard 4d geodesic equations 
(|Maartens fc Kovamal [2OI0I ) . However, the bulk has an ef- 
fect on the motion of the test particles on the brane via the 
metric. Since the projected Weyl tensor effectively serves as 
an additional matter source, the metric is affected by these 
bulk effects, and so are the geodesic equations. This has 
to be contrasted with Kaluza-Klein theories, where matter 
travels on 5d geodesies. 

In order to obtain results which are relevant to the 
galactic dynamics, in the following we will restrict our study 
to the static and spherically symmetric metric given by 



ds 



e-^^^dt^ + e^^^^dr^ + r^dn^ 



(13) 



where dO,'^ = dO"^ + sin^ ( 

The Lagrangian C for a massive test particle traveling 
on the brane reads 

C^^{-e"f + e\-^ + r^n') , (14) 

where the dot means differentiation with respect to the affine 
parameter. 

Since the metric tensor coefficients do not explicitly de- 
pend on t and Q, the Lagrangian (|14p yields the following 
conserved quantities (generalized momenta): 
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t = E, 



(15) 



where E is related to the total energy of the particle and L 
to the total angular momentum. With the use of conserved 
quantities we obtain from Eq. (|14l) the geodesic equation for 
massive particles (for which 2jC — —1 holds) in the form 



e r + e 



(16) 



The second term of the left-hand side can, in some cases, 
be interpreted as an effective potential. For instance, for 
the Schwarzschild space-time, where e"^^^ = 1, the kinetic 
term is position independent. In that case the notion of an 
effective potential is appropriate. In other cases, even one 
can still compute the turning points of the kinetic term, 
however, the effective potential interpretation is lost. 

For particles in circular and stable orbits the following 
conditions must be satisfied: a) f = (circular motion) b) 
dVeff/dr = (extreme motion) and c) d^Vef f /dr ^ \^^^^ > 
(stable orbit), respectively. Conditions a) and b) immedi- 
ately give the conserved quantities as 



E' 



and 



1 + — 



2 



(17) 



(18) 



respectively. Equivalently, these two equations can be rewrit- 
ten as 

r2 r^^'l2 



1 - rv'/2' 



(19) 



We define the tangential velocity vtg of a test par- 
ticle on the brane, as measured in terms of the proper 
ti me, that is, by a n observer located at the given point, 
as l|Landau fc Lifs hitz 1975)) 



2 -v 2A2;-2 -v^ ;-2 

Vtg = e r \l t = e —t 



(20) 



In the second equality we have employed the second Eq. (|15p . 
By eliminating L with Eq. 1)18^ and subsequently E with 
the first Eq. (|15|) . we obtain the expression of the tangential 
velocity of a test particle in a stable circular orbit on the 
brane as (|Matos et al.ll2000l : iNucamendi et al.ll20oil ) 



2 
Vtg 



TV 



(21) 



Let us emphasize again that the function v' is obtained 
by solving the field equations containing the bulk effects as 
additional matter terms; we consider this in Section [S] 



2.3 



The gravitational field equations for a static 
spherically symmetric brane 



For the metric given by Eq. (|13[) the gravitational field equa- 
tions and the effective energy-momentum tensor conser va- 
tion equation in the vacu um take the form (^Harko fc MakI 
l2003l : lMak fc Harkdl20o3 ) 



(22) 



([/ + 2P) , 



-a1 
^ 2 



,'2 



1^- v'-\' 
2 r 



U' + 2P' 
' 2U + P 



6P 



(23) 



(24) 



(25) 



ri2U + P)' 

where ' — d/dr, and we have denoted at — kf/S. Note 
that Eq. ((25l) is a consequence of Eqs. ([H]), and ((24)) . 
respectively. 

As for the motion of the test particle on the brane we as- 
sume that they follow stable circular orbits, with tangential 
velocities given by Eq. (|2ip . Thus, the rotational velocity of 
the test body is determined by the metric coefficient exp (f) 
only. 

The field equations (l22)| - ((23)) yield the following effec- 
tive energy density p"^, radial pressure P'^^ and orthogonal 
pressure Pj^ , respectively, 

= 3abU, 
P"^ = at{U + 2P) , 
Pf = aUU-P), 



(26) 
(27) 
(28) 



which obey p° 



2P'f 



0. This is expected for 



the 'radiation' like source, the projection of the bulk Weyl 
tensor, which is trace-less, = 0. 



3 STRUCTURE EQUATIONS OF THE 
VACUUM IN THE BRANE WORLD 
MODELS 

Eq. (I22|l can immediately be integrated to give 



^-A _ ^ Ct GMu (r) 



(29) 



where Cb is an arbitrary constant of integration, and we 
denoted 



GMu (r) = Scib / U{rydr. 



(30) 



The function Mu is the gravitational mass correspond- 
ing to the dark radiation term (the dark mass). For U — 
the metric coefficient given by Eq. ((29} must tend to the 
standard general relativistic Schwarzschild metric coeffi- 
cient, which gives Cb ~ 2GM, where M = constant is the 
baryonic (usual) mass of the gravitating system. 

By substituting i^' given by Eq. (I25|l into Eq. (|23p and 
with the use of Eq. (|29p we obtain the following system of 
differential equations satisfied by the dark radiation term 
U , the dark pressure P and the dark mass Mu, describing 
the vacuum gravitational field, exterior to a m assive body, 
in the brane world model ((Harko fc Mak[ [20031 '): 



dU _ 2v% {2U + P) 
dr r 



ydP__eP_ 

dr r ' 



dMi 



u 



dr 



Sab 



with the tangential velocity given as 

1 2GM + GMu + ab{U + 2P) r^ 



2 
Vtg 



, ^-^ _ 2GM 



GMu 



(31) 



(32) 



(33) 
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In order to close the system a supplementary functional 
relation between one of the unknowns U, P, Mu and vtg 
is needed. Once this relation is known, Eqs. (|3ip - p3p give 
a full description of the geometrical properties and of the 
motion of the particles on the brane. 

The system of equations (jSip and H32p can be trans- 
formed to an autonomous system of differential equations 
by means of the transformations 



acceleration expansion shear vorticity 



= In r, g = 1 



2GM GMu 



(34) 
(35) 



We shall call n and p the "reduced" dark radiation and 
pressure, respectively. 

With the use of the new variables given by Eqs. (|34p . 
Eqs. ((321) and (|3l]) become 

§ = f*-g, (36) 

S^-i^^i±^Ui±i(fi±M„4 + 2,-2p. (37) 
do 1 — ? do 

Eqs. ((ST)) and (|32j, or, equivalently, ((36} and ^7} may 
be called the structure equations of the vacuum on the 
brane. In order to close this system an "equation of state", 
relating the reduced dark radiation and the dark pressure 
terms is needed. Generally, this equation of state is given in 
the form P = P{U). 

In the new variables the tangential velocity of a particle 
in a stable circular orbit on the brane is given by 

2 1 g + Km + 2p) 

'"ts^o — , — : ■ (38) 



By using the expression of the tangential velocity, Eq. (|37 
can be rewritten as 



^ + 2p) = -2 (2^ + p) vl + 2m - 2p. 



(39) 



Eqs. (|38|l and (|39p allow the easy check of the physical 
consistency of some simple equations of state for the dark 
pressure. The equation of state /x + 2p = immediately gives 
Vtg = 1 and q = 2/3, respectively, implying that all test 
particles in stable circular motion on the brane move at the 
speed of light. This result contradicts the assumption that 
the test particles are time-like, as well as the observations 
on galactic scale. Therefore the equation of state -f 2p = 
is not consistent with the rotation curves. The equation of 
state 2fj, + p = gives fj, — fJ-o/r^, where fio = constant is an 
arbitrary integration constant, U — ^o/Sofbr^ and GMu = 
—jj,o/3r^, respectively. In the limit of large r, the tangential 
velocity Vtg tends to zero, vtg — > 0. Therefore, this model 
seems also to be ruled out by observations. The case i-i = p 
gives ^{0) = /ioexp[— 2 J Vtg(6)d6^, /io = constant, and 



(2«?,-/i)/(l + 2«|,). 



4 TANGENTIAL VELOCITY FOR GALAXIES 
FOR A LINEAR EQUATION OF STATE OF 
THE WEYL FLUID 

4.1 Linear equation of state for the Weyl fluid 



A 











Table 1. Kinomatical 


scalars of the 


congruence in 


decomposition of VqU;, 








acceleration 


expansion 


shear 


vorticity 





e 









In order to close the system of equations Eqs. (|36p and (|37p . 
we need to specify the equation of state relating the dark 



Table 2. Kinematical scalars of the congruence in the 2-1-1 
decomposition of DaTi,. 



pressure to the dark radiation. In the full 5-dimensional ap- 
proach of IWisemanI (|2002| ) the Einstein field equations were 
solved numerically for static, spherically symmetric matter 
localized on the brane, yielding regular geometries in a bulk 
with axial symmetry. For this a density profile, taken as a 
deformed top hat function, was imposed. 

An alternative approach for closing the system of field 
equations can be obtained from the 3-1-1-1-1 covariant ap- 
proach in brane worlds, developed in iKeresztes fc Gergelvl 
(|2010al lbh. First, we impose: i) cosmological vacuum in 5- 
dimensional spacetime, ii) the brane embedding is symmet- 
rical, iii) fine-tuning on the brane (the brane cosmologi- 
cal constant vanishes), iv) no matter on the brane, v) the 
brane is static and spherical symmetric, with metric given 
by Eq. (|13p . The set of gravito-electro-magnetic quantities 
in the 3-1-1-1-1 covariant approach is presented in Appendix 
I. Due to assumptions ii) and iv) the Lanczos equations 
impose that the tensorial and vectorial projections of the 
extrinsic brane curvature along the brane normal n"* van- 
ish: (Jab = = Ka- Therefore Sab = Eab — kiPab/2 and 
Hab = Hab, respectively. 

From the above assumptions we also find that the only 
nonzero kinematical quantity related to the normal it" of 
the t = const, hypersurfaces is its acceleration Aa — u^DbUa 
(here Da is the covariant derivative on the brane) . The kine- 
matical quantities of the vector congruences it" and are 
enlisted in Tables [T] and [21 

Due to the spherical symmetry, assumption v) , it is con- 
venient to decompose the brane spacetime into a 2+1-1-1 
form. Then the gravito-electro-magnetic quantities appear- 
ing in the brane equations further reduce to Hab = 0, 
Eab = E{r) {ran - hab/3), kjU {r), Qa = and kjPab = 
kiP (r) (farb — hab/S), respectively, while Aa — A{r)ra. 
Therefore under the assumption of spherical symmetry the 
Weyl fluid is characterized by the set U (r) , P (r); the elec- 
tric part of 4-dimensional Weyl tensor by E (r); and the 
acceleration of the time-like normal to the 3-space by A (r). 

We also introduce the covariant derivative Da asso- 
ciated with the 3-metric hab and the expansion Q (r) of the 
radial geodesies in the local three-dimensional space (obey- 
ing '^^^DaTb = 0) by the relation '■^^D^Va = 6 (r). The 4 
independent field equations for the variables U, P, E, A, 
and Q, arising in the 2-1-1-1-1 brane formalism are given in 
Appendix II. There are two algebraic equations [Eqs. (jlOSp 
and (|112p ] determining U and P as function of E,A, and O: 
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ktu 



4£ 
~3~ 



ktP = e 



(40) 



(41) 



and two first order non-linear ordinary differential equations 
[Eqs. (|114|l and (|115|) ] for the variables A, 9 (also containing 

Ey- 



re-, 



4E 



(42) 



(43) 



-e' + - + -+Ae = o, 

2 4 3 r2 

The last two equations contain only quantities defined and 
determined by the brane dynamics. In order to obtain a spe- 
cial solution of the system (|42|) - (|43| ). we need a third relation 
between the brane kinematic quantities A, and electric 
Weyl brane curvature E. In order to establish this, first we 
remark, that in the Schwarzschild case these quantities obey 



2E 



+ 9^1 = e 



+ A 



1 



0. 



(44) 



By virtue of these the two Eqs. H42p - (|43p are found to co- 
incide, thus A, Q and E are still defined by a set of three 
equations. For a spherically symmetic solution on the brane 
we allow for a slightly modified identity as compared to the 
first Eq. (gH): 



2E ~, ,;^/e ,\ B 

— -i-e^ = ^e - +A] - — ., 



(45) 



with A and B two constants, both reducing to 1 in 
Schwarzschild case. By employing Eqs. (|40|l . pif) . this equa- 
tion can be rewritten as a simple equation of state for the 
Weyl fluid: 

^-(«-2)t/-^. (46) 

Here we have introduced the new constants a, B by redefin- 
ing A = a/ {2a - 3), B = {a - B) / (2a - 3). We also remark 
that in the variables given by Eq. (|34|l the equation (|46|) 
takes the simple linear form 

p(M) = (a-2)M-B. (47) 

The system H42p - (|43p and (|45p . being a first order ordinary 
system of differential equations, determines the variables A, 
G and E. Their existence is assured by the Cauchy-Peano 
theorem. When rewritten in metric variables, the solution 
valid in the region where rotation curve data is available 
will be obtained in the following subsection. 

One can be rightously worried about the compatibility 
of Eq. (|46|l with the full 5-dimensional gravitational dynam- 
ics, however this is but a boundary condition imposed on 
the brane at some arbitrary time. Then the static character 
of the problem assures that it is preserved by the dynamics 
throughout the temporal evolution. 

Looking to this problem from a mathematical point 
of view, for the system of first order differential equa- 
tions describing the the extra dimensional evolution 
iKeresztes fc Gergelvl l|2010ah . the choice of the Weyl fluid 
variables represents an "initial condition" in any suitable off- 
brane parameter (for example, the one associated with the 



integral curves of the vector field n^). Such "initial values" 
could be chosen arbitrarily except when there is a constraint 
to be obeyed on the "initial surface". This surface being the 
brane, the condition to be obeyed is the constraint equation 
for the Weyl fluid, Eq. (|11|). As shown in Appendix II, on 
a static and spherically symmetric brane Eq. (|11|) simplifies 
to Eq. (|115p . It is not very difficult to see, that this equa- 
tion follows from the system (|40p - (|43p . Therefore the con- 
straint for the Weyl fluid is trivially satisfled in the static 
and spherically symetric case and the particular choice of 
"initial conditions" (|45p is allowed, as would be any other 
choice. The advantage of chosing Eq. (|45p is, however, that 
it leads to the simple linear equation of state in the reduced 
Weyl variables p and /j.. We will see that the confrontation 
of the model with galactic rotation curve data supports the 
choice (|45p from a physical point of view either. 



4.2 The metric and tangential velocity on the 
brane 

The dark radiation and the dark pressure can be obtained as 
a functions of the tangential velocity in a closed analytical 
form for the equation of state given by Eq. I|47p . The reduced 
dark radiation can be obtained as 



Q2{3-a)/(2a-3) 



exp 



2a 



2a- 



3B 



exp 



2a- 
2a 



, (6») de 



-2(3-a)/(2a-3) 



2a- 



(e) de 



(48) 



where Co is an arbitrary integration constant. Hence, if the 
velocity proflle of a test particle in stable circular motion is 
known, one can obtain all the relevant physical parameters 
for a static spherically symmetric system on the brane. 

By using the linear equation of state of the dark pressure 
Eq. (|37)l takes the form 



(2a - 3) 



dfi 



{afi - B)[q+ (2a - 3) m/3 - 2B/3] 



2(3-a)^-h2B, 



(49) 



where we have neglected the possible effect of the cosmo- 
logical constant on the structure of the cluster. Due to the 
mathematical structure of Eq. (|49p there are two cases that 
can be considered separately, a = 3/2, and a 7^ 3/2, respec- 
tively. 

For a galactic dark matter halo with mass of the order 
of M = W^^Mq and radius R = 100 kpc, the quantity 
2GM/R is of the order of 9.6 x 10~^, which is much smaller 
than one. Since observations show that inside the galaxy 
the mass is a linearly increasing function of the radius r, the 
value of this ratio is roughly the same at all points in the 
galaxy. Therefore from its definition it follows that generally 
q « 1, and 1 — q ~ 1. Moreover, the quantities q^ and 
qdq/dO are also very small as compared to q. Eq. (|36p gives 
^l^q + dq/de, dfi/dO = dq/dO + d^q/dO'^. 

4.2.1 The case a^i/2 

Hence, by neglecting the second order terms and assuming 
a 7^ 3/2, we obtain for q the following differential equation: 
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(50) 



cPq dq 
_+m--nq = b. 



a/ 



2' 



(51) 
(52) 

(53) 



where we have denoted 

^ B 2 a (B - 3) + 9 

m — 1 ^ 

3 3 2a - 3 

2 a (2B - 3) + 9 ,3 

3 2a - 3 ' ^ 2 ' 

and 

2 5(5-3) 3 
''=3^^2^'"^2' 

respectively. The general sohition of Eq. H50p is given by 
q{e) = vo + Cie'''' + C2e''',a^^, (54) 

where Ci and C2 are arbitrary constants of integration, and 
we denoted 

_ _b _ B{B-3) 
"° " "ri " a(2B-3) + 9' 

and 



(55) 



-m ± Vr 



4?! 



(56) 



respectively. The reduced dark radiation term is given by 



fi (9) = vo + Ci{l + h) e'" + C2 (1 + h) e- 



(57) 



The tangential velocity of a test particle in the "dark 
matter" dominated region is given by 



2 



or 



q+-(p + 2p) 



q + 



(2a - 3) ^ - 2B 



Vtg (d) ~ Vtgao + 76 ' + ^ , 

where we have denoted 

Vtgoo = ^ (ffl'^O - B) , 



'^^2 



2a 



3 

2a - 3 



(l + h) + 1 



(l + /2) + l 



Ci 



C2. 



(58) 
(59) 

(60) 
(61) 
(62) 



In the initial radial coordinate r the tangential velocity 
can be expressed as 



vtg (r) ~ Vtgoo + + W" 



(63) 



In order to have an asymptotically constant Vtg (r) , both 
li and I2 should be negative numbers, h < and Z2 < 
, respectively. This can be achieved if m > 0, n < hold 
simultaneously. 

In the original radial variable r we obtain for the dark 
radiation and the mass distribution inside the cluster the 
expressions 



'^~' + C2(l + /2)r (64) 



3atU{r) ^^+Ci{l + h)r 
and 

GMu{r) = r{vo + Cir' + Car'^) - 2GM, (65) 

respectively, where we have neglected the possible effect of 
the cosmological constant. The first term in the mass profile 



of the dark mass given by Eq. (|65p is linearly increasing with 
r, thus having a similar behavior to the dark matter outside 
galaxies. The second and third terms in Eq. (|65|l are either 
constants (for Zi,2 = —1) or they decrease with increasing r, 
as /i,2 < 0. 

Finally, the metric coefficient f can be calculated from 
the equation 



du{e) 2a -3 2B 



(66) 



givmg 



Cr^^'s- X 



exp 



■ 3 + (2a-3)(l + /i) ^,, 

3h + 



^ 3+(2a-3)(l + /2) 
3h ' 



(67) 



where is an arbitrary integration constant. In the limit 
of large distances e"^^^ behaves like e"'""' « C^r'^"'!'°° . Inside 
the galaxy we can approximate ~ 1 — C/r — GMu{r) jr. 

A few comments are in order. First, in the particular 
case i? = 0, we have uo = 0, and in consequence the obtained 
mass profile is not consistent in general with observations, 
as for such an equation of state the dark radiation usually 
cannot play the role of the dark matter (it has no r~ term). 
For example in the case a = 4, corresponding to an equation 
of state P = 2(7, the dark mass is given by GMu(r) — 



Cb. 



Secondly, the effective, geometry induced mass must 
satisfy the condition Mu ^ 0. Therefore the solution ob- 
tained in the present Section is valid only for values of the 
coordinate radius r so that vor + Cir'^"'"^ + C2r'2+^ — G^O. 
In the limit of small r , taking into account that Zi < and 
I2 < 0, and by assuming that the coefficient of the dominant 
term (either vo, Gi or G2) is positive, the dark mass diverges 
at the center of the cluster, limr^o Mu{r) = 00. Otherwise 
there is a point ro, where Mu{ro) ~ and Muir) < 0, for 
r < r-Q. In this case the solution has physically acceptable 
properties only in the region r ^ ro. 



5 FURTHER SIMPLIFICATION FROM 
OBSERVATIONAL CONSTRAINTS 

In Section |4] the assumption q was made in order to de- 
rive the analytical solution. As it can be seen from its defini- 
tion, he quantity q is basically a post-Newtonian parameter 
and as such the assumption is justified when studying galac- 
tic rotation curves. Further, as we are studying bounded mo- 
tions in the galaxy, due to the virial theorem should be 



of the same order. Eq. (|38j, and the remark preceding Sub- 
section IV. A , stating <C 1, together imply p <C 1. Then, 
from the definition of p, Eq. (|47[) the condition B <^ a — 2 
stems out. In other words, B/a should be of the same order 
of magnitude as /i and q, respectively. 

In what follows, we will use the above remark for sim- 
plifying the generic results established for a 7^ 3/2 in the 
previous Section. For the parameters m, n we obtain 

4a - 9 „ a - 3 



III' ~ ~ , 

2a - 3 ' 
which give 



2a - 3' 



(68) 
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h 



-1 + 



(69) 



2a - 3 

In order Z2 < to hold (or equivalently m > and n < 0), 
the constant a must take values in the range a G (—00, 3/2)U 
(3,00). 

With these values of Zi,2 the constants C1.2 simplify to 
7 Ci/2 and « C2. We also get 



vo 



B{B~3) /a 



B 



(70) 



2B-3 + 9/a 

for any a in the allowed domain. The rotational velocity at 
infinity then becomes 

a / B\ B 
3 



[vo ~ — 

V a / a - 



(71) 



Obviously Vtg^ ^ 0, therefore either a < 3/2 and B ^ 
or a > 3 and B > 0. 

The metric functions approximate as 



Q rr^ \ ^ ""0 , 3C2 -3(1- 



GMu{r) « vor + Ci + C2r2o- 



2GM, 



G^r '9°= exp 



-1 „ 2a — 3 
dr —02 —r ' 



(72) 

(73) 
(74) 



a — 3 

and for a generic r the rotational velocity can be written as 



vig {r) 



r + — r +C2r ^^"-3. 

a — 3 2 



(75) 



By defining a = 3/ (2a - 3) = Z2 + I and /3 = B/ (a - 3), 
respectively; then re-introducing the velocity of light c for 
dimensional reasons; by replacing Ci with a mass type con- 
stant M-\-Mq= C\/2Gc? , and also by rescaling the constant 



C2 as C2 = C(?rl " , we obtain 



+I3+C 



, r > r*.(76) 



Here we have changed the notations (M, Mo) to 
(Mj°',Mo°*) for reasons to be explained in what follows. 
This solution is valid for any r > r*, where r* represents 
the radius beyond which the baryonic matter does not ex- 
tend. We interpret the constant Mj°* as the total baryonic 
mass inside radius r* , while Afo°* is a universal constant 
due to the brane fluid. Thus the Weyl fluid is characterized 
by three universal dimensionless constants a, (3, C and a 
mass- type constant Mo°*. Also an arbitrary scaling constant 
rt, was introduced in order to have C dimensionless. 

Due to the constraints established for the Weyl fluid 
parameters a and B , the newly introduced constants a, /? 
obey either a<QorQ<a<l (the value a = is excluded, 
as it would correspond to the unphysical values a — >■ ±00) 
m both cases P being a small positive number < /? ^ 1 
(with the exception of a ~ 3, translating to a ~ 1, when j3 
can be an arbitrary positive number). 



6 HIGH SURFACE BRIGHTNESS GALAXIES 

6.1 The baryonic sector: the bulge-disk 
decomposition 

We model the distribution of baryonic mass in High Surface 
Brightness galaxies as a sum of disk and bulge components 



with constant, but distinct mass-to-light ratios. We estimate 
the bulge parameters from a Sersic r^^" bulge model and the 
disk parameters from an exponential disk model, both fitted 
to the optical I-band galaxy light profiles. 

The surface brightness (specific intensity) of the 
spheroidal bulge componen t of each g alaxy is given by a 
generalized Sersic function l|Sersidl 19681 ): 



h{r) = /o,6exp 



l/n 



(77) 



where Jo, 6 is the central surface brightness of the bulge, ro is 
its characteristic radius and n is the shape parameter of the 
magnitude-radius curve. As a rule, early-type spiral galaxy 
bulges have n > 1 , while late-type spiral galaxy bulges are 
characterized by n < 1. The radius of the bulge Vb is defined 
by the condition of the surface brightness being equal to 
^0 = 2, 64 X 10"* mjy- arcsec"'^ in the I-band images. 

In a spiral galaxy, the radial surface brightness profile of 
the disk exponentially decreases with the radius (|Freemanl 
1970) 



(78) 



Idir) = /o,d exp ( 



where 7o,d is the disk central surface brightness and h is 
a characteristic disk length scale. In order to measure the 
light and mass distribution in the disk, the model image of 
the bulge is subtracted from the original I-band image. All 
remaining light in these images is then assumed to originate 
from the disk component. 

The (bolometric) luminosity of any of the galaxy com- 
ponents is an integral over surface, solid angle and fre- 
quency of lb and Id, respectively. The respective mass over 
luminosity is the mass-to-light ratio - for the Sun 79 = 
5133 kg W~^. The mass-to-light ratios of the bulge and 
disk a and Tb will be given in units of 7© (solar units). 
We will also give the masses in units of the solar mass 
Mq = 1.98892 X 10^° kg. Assuming that the mass dis- 
tribution of a spiral galaxy follows the de-projected sur - 
face brightness distribution dK annapan fc Gawiseij |20o3); 
IPortinari et al. (12004) ) (corrected wit h respect to the incli- 
nation, cf. (jPalunas fc WilUamslbOOOl ). with constant mass- 
to-light ratios a and Tb the bulge and disk masses within the 
radius r are 



Mb{r) = a^^^2Tr I Ib{r)rdr, 
Iq 



Md{r) = rb^^^P-27v I Id<ir)rdr, 



(79) 
(80) 



where Fq (D) is the apparent flux density of the Sun at a 
distance D Mpc, Fq (D) = 2.635 x 10'^-°"-^Q mJy , with 
/q = 4.08 + 5lg {D/1 Mpc) + 25 mag, and 



Af{D) = 4.4684 x 10~^^i3"^ m" 



(81) 



The baryonic matter at any r is given by 

Mbaryonic{r) = Mb{r) + Md{r). (82) 

In iPalunas fc Williams! l|2000l ) a maximum disk mass 
model was presented, and the bulge-disk decompositions was 
performed for a sample of 74 high surface brightness spiral 
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galaxies, for which the I-band surface brightness profiles and 
Ha velocity profiles were also given. From this set we have 
left out those galaxies, which either have bars or rings, which 
would both contradict the assumption of spherical symme- 
try. In addition we have left out the galaxies for which the 
above baryonic model can not be applied, since they have 
non- vanishi ng bulge surface brightness values, but vanishing 
bulge mass (|Palunas fc Williamsl '2000Y And finally we have 
also omitted the galaxies where the observed data sequence 
exhibits a wavy pattern instead of a plateau, suggesting a 
disk structure strongly contradicting our assumption for ro- 
tational symmetry about the galaxy center. Such patterns 
are illustrated in Fig. [1] As a result of this selection pro- 
cess we end up with 9 galaxies. Their rotation curves are 
represented on Fig. O 

As for the bulge parameters of the 9 selected galaxies, 
we have derived the best fitting values of 7o,i,, n, ro, rt, 
7o,d, and h from the photometric data as well as the mass- 
to-light ratio of each component a and rj, by fitting Eq. 
II82I) to the data on rotatio n curves represented on Fig. 1 of 
iPalunas fc Wilhamj (|2000D . These are all collected in Ta- 
ble O 



6.2 Combined baryonic and Weyl model 

We assume that within the bulge radius rj, the contribu- 
tion of the Weyl-fluid can be neglected and the part of the 
observed rotation curves lying below rt could be explained 
with baryonic matter alone. Outside this radius we switch 
on the Weyl-fluid and we take into account the exponential 
disk as a perturbation that does not affect the geometry. We 
add, however, its contribution to the rotation velocity. 



,M = 



G[Mtir) + AUr)] 



+ 



P + C 



Hk in) , 



(83) 



where Hk (rt) is some function smoothly approaching the 
Heaviside step function: 



H{n)^ lim Hk{n) = 

k >-oo 



r <rb 
r ^ rt 



(84) 



The parameter Mq"* of the Weyl fluid in Eq. ((76} repre- 
sents the contribution of the equivalent Weyl mass from the 
spheres within a given radius. As such a contribution is taken 
care by the first terms of Eq. (|83|l . we have chosen Mq°* = 0. 
For Hk (rt) we choose the logistic function: 

^-(")- l + exp(-L(r -..)) • ^''^ 
The value of k gives the sharpness of the transition. We 
determine the parameter k by imposing the conditions 
Hk (0.95ri,) = 0.001, which also implies Hk (LOSr;,) = 0.999 
and Hk [rt) = 0.5 and its val ues are given in Table [5] for 
each of the 9 selected galaxies ()Palunas fc Williamsll2000i l. 

The baryonic parameters r;,, a and the Weyl parame- 
ter a have to be determined by a ^^-minimization fitting 
of Eqs. (|77p - H85p . with the rotation curve data for each 
individual galaxy. The parameter /3 of Weyl sector gives the 
asymptotic rotation velocity. We identify /3 with the average 
of the points of the rotation curves for r > r% (these are on 
the plateau). 



At the end of this subsection we prove that the "trun- 
cation" of the Weyl fluid at rb does not induce any distribu- 
tional source layer at rj, thus it is consistent with the junc- 
tion conditions across the sphere with radius rb- For this we 
first remark that at a formal level the truncation can be im- 
posed in all equations by the replacements /3 — >■ fiHk (rb). By 
also keeping in mind that C = C2C~^r'^~^ — —fi was chosen, 
Eq. (|72p shows that U smooths to zero with Hk (vb) across 
n. The equation of state (gZ} and B = (3/3/2) (1/a - 1) 
then guarantees that P has the same property. As both in- 
side and outside rb we have the same spherically symmetric 
metric and the energy-momentum tensor changes smoothly 
across rb, there is no discontinuity in either the metric or 
its first derivative, therefore the even stronger Lichnerow- 
icz continuity conditions are obeyed. The Israel-Lanczos- 
Darmois junction conditions of the continuity of induced 
metric and extrinsic curva ture thus f ollow . Our approach is 
different from that used bv lWisemanI (|2002l l. in which a given 
density profile is considered, together with the condition 
of the isotropy on the brane. Therefore, while in IWisemanI 
([2OO2I) the matching conditions impose these requirements 
as boundary conditions, in our approach the matching con- 
ditions are automatically satisfied. 



6.3 HSB galaxy rotation curves 

Despite of the differences in the surface brightness profiles 
and rotation curves for the chosen HSB galaxies, the com- 
bined baryonic-|~Weyl model Eq. (|83p fits well the sample. 

A restriction on the parameter space emerges as 13 — 
— C from the condition that the Weyl contribution to Vtg (r) 
vanishes at rt, (by switching on the Weyl fluid at r = Vb). 
Therefore Eq. (|83|l simplifies to 



vtg{r) 



G[Mb{r) + M4r)] 



2 o 

c /3 



(7)" 



+ 



Hk (rb) 



(86) 



The rotation curves themselves together with the data are 
shown on Fig. [2l 



7 LOW SURFACE BRIGHTNESS GALAXIES 
7.1 LSB galaxy rotation curves 

A typical LSB galaxy resembles a normal late-type spiral, 
usually with some ill-defined spiral arms. They usually have 
HI masses of a few times 10^ . iThuan et al.l l|l987l ) and 
iBothun et all (| 19931 ') have found that, although LSB galax- 
ies follow the spatial distribution of HSB galaxies, they tend 
to be more isolated from their nearest neighbors than HSB 
galaxies. LSB galaxies can be distinguished from the galax- 
ies defining the Hubble sequence by their low surface bright- 
ness, rather than small size (LSB galaxies are not necessar- 
ily dwarfs). The central surface brightness of LSB galaxies 
is much lower than /ifl(O) — 21.65 ± 0.3 mag-arcsec^ - the 
typical B-ban d value for HSB galaxies, established by the 
Freeman law jFreemanl l|l970l ): ?). 

The mass density distribution of LSB galaxies at small 
radii is dominated by a nearly constant density core with 
a total mass Mq, and a radius rc of only a few kpc, as 
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ES0215G39 


61.29 


0.1171 


0.6609 


0.78 


2.58 


0.0339 


4.11 


26.28 


0.04 


2.1 


0.67 


2.37x10 


-7 


28.29 


ES0322G76 


64.28 


0.2383 


0.8344 


0.91 


4.50 


0.0251 


5.28 


15.35 


0.47 


3.28 


0.76 


3.21x10 


-7 


38.69 
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38.19 


0.1949 


0.7552 


0.33 


1.37 


0.0744 


2.20 


50.32 


1.27 


2.7 


0.57 


3.92x10 


-7 


10.15 
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0.1113 


0.4626 


0.43 


0.99 
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3.47 


34.58 


2.96 


1.85 


0.69 


5.25x10 


-7 


34.77 


ESO383G02 


85.40 


0.6479 


0.7408 


0.42 


1.94 


0.5118 


3.82 


17.79 


0.47 


0.22 


0.70 


3.72x10 


-7 


21.55 


ES0445G19 


66.05 


0.1702 


0.6133 


0.57 


1.79 


0.0478 


4.27 


38.59 
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1.39 


0.46 


3.64x10 


-7 


29.26 


ESO446G01 


98.34 


0.2093 


0.8427 


1.28 


6.33 


0.0357 


5.25 


10.90 


0.93 


2.82 


0.54 


4.64x10 


-7 


43.35 


ESO509G80 


92.86 


0.2090 


0.7621 


1.10 


4.69 


0.0176 


11.03 


14.75 


0.61 


5.50 


0.87 


6.3x10- 


7 


25.71 


ES0569G17 


57.77 


0.2452 


0.4985 


0.45 


1.18 


0.1348 


2.06 


58.74 


0.06 


1.4 


0.65 


3.18x10 


-7 


7.24 



Table 3. The baryonic parameters (D, /o,iji '^i ''0, ''iji ^odi h, k) of the 9 HSB galaxy sample. The additional baryonic parameters cr, rj, 
and the Weyl parameter a are determined by fitting. The best fit parameters (cr = Mi,/Li,, ti, = Mj^/L^, a) and the minimum value 
of the statistic Xmin ^^^'^ given. All Xmin ^"'-^ within la. 
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DDO 189 


57.5 


4.05x10** 


1.25 


0.3 


6.43x10' 


-8 
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NGC 2366 


46.0 


1.05x10^ 


1.47 


0.8 


1.12x10' 


-7 


2.538 


NGC 3274 


138.1 


4.38x10* 


0.69 


-0.4 


6.73x10- 


-8 


18.099 


NGC 4395 


30.0 


2.37x10* 


0.71 


0.9 


3.43x10" 


-7 


27.98 


NGC 4455 


99.7 


2.26x10* 


1.03 


0.9 


2.72x10" 


-7 


7.129 


NGC 5023 


86.3 


2.69x10* 


0.74 


0.9 


4.53x10" 


-7 


10.614 


UGC 10310 


36.4 


1.28x10^ 


2.6 


0.4 


1.12x10" 


-7 


0.729 


UGC 1230 


15.3 


3.87x109 


3.22 


-1.7 


1.12x10" 


-7 


0.539 


UGC 3137 


34.5 


5.32x10^ 


3.87 


-0.5 


1.23x10" 


-7 


4.877 



Table 4. The best fit parameters of the 9 LSB galaxy sample (Mq, rc , a, p). The minimum values of the statistic Xmin ^'^^ 
given. All x^ • are within la. 
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Figure 1. Three examples of galaxies omitted due to the wavy pattern data sequence. The best fits of our model are also shown, all 
outside the la confidence level. 
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ES0215G39 


1.27 


1.89 


ES0322G76 


2.43 


8.53 


ES0322G77 


2.19 


1.487 


ES0323G25 


3.75 


3.76 


ESO383G02 


2.35 


3.12 


ES0445G19 


2.33 


1.97 


ESO446G01 


4.69 


9.7 


ESO509G80 


10.7 


205 


ES0569G17 


1.55 


1.17 



Table 5. The calculated halo masses for the investigated sample of HSB galaxies. 



12 Gergely, Harko, Dwornik, Kupi, Keresztes 




Figure 2. Best fit curves for the chosen HSB galaxies compatible with the spherical symmetry assumption, the baryonic model and 
having sufliciently accurate rotation curve and photometric data available. In the baryonic matter dominated radial range r < rt, the 
model curves bend as determined by the photometric data, shown in Table [S] together with the Weyl parameters of the model. The 
discontinuity in the first derivative of the model curves indicate the bulge radius ri,. 



established by Ide Blok et all (|200ll ). Hence we ignore the 
mass contributions of the stellar and gas components. At low 
radii r < rc, the constant density implies w^(r) = GpV/r, 
where V is the volume of a sphere with radius r, and 
p = 3Mo/4rc7r, i.e., n ~ r. The Weyl fluid can reproduce 
this behavior if we choose the equation of state pip) ~ 
(i.e. a — 2 and B = in Eq. H47p . Then the dark anisotropic 
stress/pressure also vanishes (P — 0), and the remaining 
part of the Weyl fluid energy momentum tensor (as = 0) 
then represents radiation. Using Eqs. (|51l) - (|52l) . (I55p - (|56p 
and (|60p - (|62p . for C2 ~ the tangential velocity, given 
by Eq. (|63p . yields to Vtg ~ jr^ , the desired behavior. For 
r > Tc we use Eq. (|76|) . without baryonic matter (Mj°' — 0), 
and apply a new notation: Afo°* — Ma. Requiring the conti- 
nuity of the rotation velocity through r = rc we have 



Hk {',) . (87) 



where Hk (rc) is the logistic function given by Eq. (|85p . 

We test the model with a sam ple of 9 LSB galax i es, ex - 
tracted from a larger sample from Ide Blok fc Bosmal (|2002l l 
as typical galaxies exhibiting the plateau region. We fit 
Eq. (|87|) with the rotation curve data taken from combined 
HI and Ha measurements. The fitted curves are represented 
on Fig. [3l In all cases we find remarkably good agreement 
between the model and observations. Outside the core ra- 
dius the parameters of the equation of state a and B are 
determined by fitting through a and /3 (see Section 5). 



From a x^-test we determine Mo and rc, and also the 
Weyl parameters a and /3, shown on Table |4] and on Fig. [l] 



8 DISCUSSIONS AND FINAL REMARKS 

The shape of the observed rotation curves strongly indicate 
the need for dark matter or equivalent modifications of grav- 
ity on galactic scale and above. However there are no uni- 
versally accepted candidates explaining the whole amount of 
dark matter needed for agrement with observations. Because 
dark matter does not interact with ordinary matter, except 
gravitationally, the question comes whether dark matter is 
rather an effect of modified gravity. We have investigated 
whether extra dimensions could provide this type of mod- 
ification in the simplest, codimension one brane-world sce- 
nario. We have derived an analytical expression for the ro- 
tational velocity of a test particle on a stable circular orbit 
in the exterior region to a galaxy, with Weyl fiuid contribu- 
tions (the dark pressure and dark radiation) included. For 
this we have assumed a linear equation of state for the Weyl 
fluid and we have argued for the rightness of this as follows. 

In subsection 4.1 we have derived a system of two first 
order differential equations for three brane variables with 
straightforward geometrical meaning. These variables are: 
(1) the acceleration of the normal congruence to the t=const 
hypersurfaces; (2) the expansion of the normal congruences 
to the r=const spheres in the t=const hypersurfaces; and (3) 
the scalar characterizing the electric part of the Weyl cur- 
vature of the brane (due to its construction, this is defined 
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Figure 3. Best fit curves for tlie LSB galaxy sample. The parameters of the model velocity curves are given in Table |4] 



on the t=const hypersurfaces) . The dark pressure and dark 
radiation do not enter these equations. 

We also derived additional algebraic relations for these 
three quantities valid in Schwarzschild space-time. As the 
brane- world spherically symmetric space-time is a modifica- 
tion of the Schwarzschild solution, we have slightly mod- 
ified one of these relations (by introducing two continu- 
ous deformation parameters, which take the value 1 for 
Schwarzschild), closing in this way the system of equations 
for these variables. The Cauchy-Peano theorem then proves 
the existence of a solution for this system. 

We have shown then, how the dark pressure and dark 
radiation are determined algebraically in terms of the above 
mentioned two kinematical and one brane Weyl quantities. 
We find that the above mentioned choice induces a linear 
relation between them, the desired equation of state. 

There are two types of evolutions, which in principle 
could destroy this solution. No problem is posed by the tem- 
poral evolution, due to the static character of the space-time. 
On the other hand, the above determined solution of the 
brane equations stands as a set of "initial" conditions for the 
sytem of first order off-brane evolution (propagation) equa- 
tions (for example along some suitably defined parameter 
associated to the brane normal). This in turn again allows 
for solutions due to the Cauchy-Peano theorem. Thus we 
only have to check then, whether such initial conditions are 
allowed, in other words whether there is any constraint in- 
volving these initial conditions, which has to be obeyed. Such 
a constraint indeed exists: it is the 4-dimensional (twice con- 
tracted) Bianchi identity on the brane. This identity however 
trivially follows from the equations solved before, therefore 
the linear equation of state is compatible with both types of 
evolution. 

The model has several free parameters, which in prin- 
ciple could be fixed in such a way to explain the observed 
galactic rotation curve behavior. 



In order to much closely test this assumption, we have 
considered a sample of 9 HSB and 9 LSB galaxies with 
well measured combined HI and Ha rotation curves. Since 
LSB galaxies ar e known to be dark matter dominated 
jMcGaugh et alj (|200ll )). reproducing their rotation curves 
without the need of any dark matter component would be 
a major achievement. Fitting the model to rotation curve 
data and photometric measurements allowed us to constrain 
the Weyl parameters a and /3 — —C; also determine the 
mass-to-light ratios AI/L of the baryonic components in 
LISB galaxies, and the total core mass Mq and radius Vc 
of LSB galaxies. The fit was in all cases within la confi- 
dence level, which supports the choice of the equation of 
state (|46p from a physical point of view. In all cases we 
found l3 <^ 1 and in most cases a G (0.4, 0.9), which for the 
constants A = {1 + a) /2, B = A {1 - P) + a/S introduced 
in Section 4.1 give A ^ B € (0.7,0.95). These express the 
difference of our model from a Schwarzschild solution, for 
which A = B = 1. 

With the parameters determined from the fit the theo- 
retical rotation curves will have an almost fiat (slightly in- 
creasing) asymptotic behavior at larger radii then the avail- 
able observations of HI and Ha velocities for these galax- 
ies. This tendency is somewhat contradictory with the shape 
of the universal rotation curve, scaling only with the virial 
mass (|Salucci et al.ll2007l ). at least for the well-fitting model 
parameters. The asymptotic shape of the universal rota- 
tion curve shows a decreasing tendency, as has been found 
from N-body simulat ions, which assume t he Navarro-Frenk- 
White CDM model (|Navarro et al.lll996l '). We note though 
that parameters of the Weyl fiuid reproducing the asymp- 
totics of the universal rotation curve can be found, but for 
these the fit with the HI and Ha velocity data falls outside 
3(T confidence level. 

Observationally, the galactic rotation curves remain fiat 
to the farthest distances that can be observed. On the other 
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hand there is a simple way to estimate an upper bound for 
the cutoff of the constancy of the tangential velocities. The 
idea is to consider the point at which the decaying density 
profile of the dark radiation associated to the galaxy be- 
comes smaller than the average energy density of the Uni- 
verse. Let the value of the coordinate radius at the point 
where the two densities are equal be i?™^". Then at this 
point 3qi,C/ (i^g"*") = (SttG/c^) puniv, where punivC^ is the 
mean energy density of the Universe. In the limit of large 
r, with the use of Eq. (I72|) . we can approximate the dark 
radiation term as 3a(,(/ ~ vo/r'^. Elence for Ru we obtain 



7->max 

Ru = c 



univ 



The mean density of the Universe is given by puniv ~ 
Pcrit = 3Hq /8tvG, where Ho is the Hubble constant, given 
by Ho = WOh km/sec Mpc, 1/2 < ft 1. Therefore 



7-) max 

Kjj 



cm 



3 X 10^ 



Mpc. 



(89) 



A numerical evaluation of R^^^ requires the knowl- 
edge of the constant vo in the flat region, and of the 
basic fundamental cosmological parameters. In the case 
of the HSB galaxies vo G (0.237 x 10"^ 0.63 x 10"^) , 
while in the case of the LSB galaxies, vo G 
(0.316 X 10-^ 0.316 X 10"'^). This gives for a range of 
Ru'^'' e (0.84/i~^ Mpc, 1.37/i"^ Mpc) for the HSB galaxies, 
and Tig"'" G (0.307/1^^ Mpc, O.QTft-^ Mpc) for the LSB 
galaxies, respectively. The measured flat regions are about 
R 2x Ropt , where Ropt is the radius encompassing 83% of 
the t o tal integrated light o f the galaxv (pinnev fc Tremaing 
1 19871 : iPersic et al.l Il996l : iBoriello fc Saluccil bOOl*). If we 
take as a typical value 7? ~ 30 kpc, then it follows that 
R « RYr^^. However, according to our model, the flat 
rotation curves region should extend far beyond the present 
measured range. 

An alternative estimation of R^^^ can be obtained from 
the observational requirement that at the cosmological level 
the energy density of the dark matter represents a fraction 
Qm ~ 0.3 of the total energy density of the Universe f2 = 1. 
Therefore the dark matter contribution inside a radius i?™^" 
IS given by 47rt7„ {Rf^^'^'f Pcr^t/S, which gives 



I 1 C ^Vtg_\ 

\l 2nmHo \ c J 



(90) 



Therefore, by assuming that the dark radiation contri- 
bution to the total energy density of the Universe is of the 
order of fi^ ^ 0.3 we have RY^""" G (388/i"\ 3881/i"^) kpc 
for vtg G (lO-^10-^). 

The limiting radius at which the effects of the extra- 
dimensions extend, far away from the baryonic matter dis- 
tribution, is given in the present model by Eqs. (|89|l or (I90p . 
In the standard dark matter models this radius is called the 
truncation parameter s, and it describes the extent of the 
dark matter halos. Values of the truncation parameter by 
we ak lensing have been o btained for several fiducial galaxies 
bv iHoekstra et all Hooi). In the following we compare our 
results with the observational values of s obtained by fitting 
the observed values with the trunc ated isothermal spher e 
model, as discussed in some detail in lHoekstra et all l|2004 ). 
The truncation parameter s is related to R^}"'^ by the re- 



lation s = RY}""^ /2-K (see Eq. (4) in lHoekstra et al.l l|2004l 
Therefore, generally s can be obtained from the relation 



2f^m 



(91) 



where a is the velocity dispersion, expressed in km/s. Hence 
the truncation parameter is a simple function of the veloc- 
ity dispersion and of the cosmological parameters only. For 
a velocity dispersion of a = 146 km/s and with Q.m = 0.3, 
Eq. (|9ip gives s ~ 245/i~^ kpc, while the t runca tion size 
obtained observationally in iHoekstra et al.l (l2004h is s = 
213/i"^ kpc. For cr = 110 km/s we obtain s « 184/i~^ kpc, 
while CT = 136 km/s gives s « 228/1" ^ kpc. All these val- 
ues are c onsistent with the observational results reported in 
:Hoekstra et al. (20041), the error between prediction and ob- 
servation being of the order of 20%. We have also to mention 
that the observational values of the truncation parameter de- 
pend on the scaling relation between the velocity dispersion 
and the fiducial l uminosity of the galax y. Two cases have 
been considered in lHoekstra et akl (|2004h . the case in which 
the luminosity Lb does not evolve with the redshift z and 
the case in which Lb scales with z as Lb <x z). Depend- 
ing on the scaling relation slightly different values of the 
velocity dispersion and truncation parameter are obtained. 

One of the most straightforward evidence for dark 
matter comes from th e radial TuUy-Fisher relation 
l|Yegorova fc Saluccil I2OO7I ): at a given galactocentric dis- 
tance (measured in unit of the optical radius) there is a 
relation between rotation velocities and the absolute mag- 
nitudes of the galaxies. From this relation we can extract 
information about the mass distribution of spiral galaxies. 
We have studied whether the Weyl model satisfies this rela- 
ti on. The plot on Fig. [Sjisjn good agreement with the result 
of lYegorova fc Saluccil (|2007l ). 

Statistical studies of N-body simulations reveal that 
there is a link between halo masses and galaxy proper- 
ties (e.g. l uminosity, baryonic mass) (|Shankar et al.l l2006l : 
iMoster et al. 2010). We have compared the equivalent halo 
mass predicted by the Weyl model Mwe.yi with the halo mass 
arising from the simulations, Msimuiation (|Moster et al.l 
l2010t ) (jSj. From the baryonic mass Mi,(oo)-|-Md(oo) we have 
derived the p redicted halo mass Msimuiation with Eq. (2) in 
iMoster et al] ||2010i ). With this mass we have calculated the 
virial radius of the galaxy Rvir (see Eq. (4) in lSalucci et al.l 
(2007*)) and the Weyl halo mass, defined as 



M, 



Weyl 



2 a 

C P 



1 



Rv 



G 



(92) 



With one exception the emerging masses are of same order 
of magnitude, in one case being also equal. These derived 
data signals some tension between the Weyl fluid model and 
the numerical simulations. 
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APPENDIX I: 3+1+1 DECOMPOSITION OF 
THE 5D WEYL CURVATURE 

With the brane normal and the temporal normal 
to the spatial hypersurfaces singled out, the five- 
dimen sional Weyl tensor Cabcd ad mits a 3+1+1 decompo- 
sition teeresztes fc Gergelvl [ioiOal ') . generalizing the corre- 
sponding decomposition of the four-dimensional Weyl tensor 

Cabcd- 

The projections 

— kiu = CABCDTL^u^ri^u^, (93) 
fclQif = CABCDhi/'n''u'^n'' , (94) 
-k^PAB = CABCohi^j^n^hi^fn'^ . (95) 

of the 5-dimensional Weyl tensor combine to an effective 
fluid on the brane, the Weyl fluid. (Angular brackets ( ) on 



indices denote tensors which are projected in all indices with 
the 3-metric hu — '^^'^ gij — ninj + uiuj, symmetrized and 
trace- free.) 

The projections 

£kl = CABCohi^j^u^h^'^u'^ (96) 

and 

T-Lkl = 2^{K ^^h]Jf Cabcdu^ (97) 

are related to the electric and magnetic parts of the four- 
dimensional Weyl tensor on the brane, Eac ~ Cabcdu'u^ and 
Hki ~ |e(fc " hj^'^CabcdU , as 




+ ^f?c(a(?6>'' + ^K(aKt) , (98) 
Hab = Hat — £(a '^'^<^h)cKd , (99) 

Here^J? = u'^u'^VcriB, Ka hJ^u^VcriK, Q = 
h^^\/ AnB, and ctab ~ hi^J^ hg'^V ktil are kinematic quan- 
tities emerging as various projections of the extrinsic curva- 
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Figure 5. The radial TuUy-Fislicr relation of our HSB and LSB 
Galaxy samples: absolute magnitude of the galaxies as function 
of the rotational velocities predicted by the Weyl model at 0.2 
(red), 0.4 (yellow), 0.6 (green), 0.8 (cyan), 1 (navy), 1.2 (lilac) 
optical radii. 



ture of the brane (0 and gab being the expansion and shear 
of the brane normal), while Va is the covariant derivative 
in the 5-dimensional space-time. 
The rest of the components 

A BC D 



£k 
Uk 

J-KL 



= CabcdIt-k u n u 

1 ab7^ c d 

— -Ek Cabcdu n , 

, A B , C D 

— CABCDn(K U l^TA n 



AB , C pi D 
£{K flL) (-'ABCDtl 



(100) 
(101) 
(102) 
(103) 



do not enter the equations on the brane. 



APPENDIX II: DYNAMICS ON THE BRANE 

In the 2+1+1 formalism on a brane with spherical symmetry 
the covariant derivative of any scalar field / (r) , taken along 
the integral c urves of the radial vector field is defined 
l|Keresztes fc Gcrgclv 2010a) as /* = r°'Daf ■ The nontrivial 
dynamical equations are four ordinary differential equations: 



e* + 



+ 



2E 



+ _± (2(7 + P) = , 



[U + 2PY +4AU + 2A + 39 P = 



A* +A e + A 



ktu = o 



yl* + ^( A-|Vi5+|p = Q 



(104) 
(105) 
(106) 

(107) 



for the five variables O, A, E, U, P . Eq. (|105|1 is equivalent 
with the constraint Eq. (fTTj) . The difference of Eqs. (|106p 
and p07|) gives a simple algebraic relation between the dark 
radiation U, dark pressure P acceleration A of temporal 
normals, expansion of radial geodesies and electric part 
of the 4-dimensional Weyl tensor E, respectively: 



kt{2U + P) = 3Ae + 2E. 



(108) 



This enables us to eliminate the Weyl constributions from 
Eq. (|106|) . obtaining 



0. 



(109) 



Another algebraic relation emerges as follows. From (a) 
the relation between the Ricci scalar of the sphere and the 
Riemann tensor of the 3-space: 



"7^ 



/ 7 ac a c\ / 1 bd b d\ 

(h — r r ) [n — r r j 



Ti-ab 



+ {Dare) (Dbra) - (Dar^) (D^r,) 



kt {U ~P)- 2E 



+ 



(110) 



where we have used DaTb = O {hat — rgTb) /2, (b) the 3- 
dimensional Riemann tensor given in iKeresztes fc Gergelvl 
(|2010al 'n. which in the spherically symmetric case simplifies 
to 



T^ab 



U heyah.b]d — 2 {^Ed[ahb\a — Ec[ahb\d^ 



^kj {Pd[ahh 



^c[ahb]d) 



(111) 



and (c) the Gaussian curvature of the two-dimensional 
spacelike group orbits orthogonal to r" and it" being TZ — 
2/r^, we obtain 



(112) 



We need to relate the newly introduced curvature coordinate 
to the ★-derivative. For this we take the ★ -derivative of Eq. 
(fTT2)) . employ Eqs. (fT09)) . ([TOSj, (fT06)) : also the ^-derivative 
of Eq. (llOSfl for eliminating _E*, finding 



Inr^ 



e. 



(113) 



This relation allows to replace the *-derivative in all equa- 
tions by r- derivatives, denoted by a prime. Thus Eqs. (|109p . 
(|105p . (|106p can be rewritten as 



4£ 



46 = 0, 



— e' + — + — 

2 2 3 

^(17 + 2P)' +4AI7 + [2A + 30^ 



re 



A' + A(@ + A 



ktU = Q 



(114) 
(115) 
(116) 



where the prime denotes the derivative with respect to r. 

In summary, the variables [/, P, A, E and Q are con- 
strained by two independent first order ordinary differential 
equations (|114p . (|116|) and two algebraic equations (|108p . 
(fTT2)) . Eq. (fTTSj) follows from these. 



^ We note that there is a missing factor 2 in front of 
S in Eq. (23) of IKeresztes fc Gergelvl l l2010al '). The conver- 
sion of notations to the notations of the present paper is: 

(£, 8a, tab, Eab) {-kjU, kjQa, - k^Pab, Eab) ■ 
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